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i Abstract 

^ I The Hartree-Fock exchange operator is an integral operator arising in the Hartree-Fock 

O^. method and replaced by a multiplicative operator (a local potential) in Kohn-Sham density 

functional theory. This article presents a detailed analysis of the mathematical properties of 
various local approximations to the nonlocal Hartree-Fock exchange operator, including the 
^ ' Slater potential, the optimized effective potential (OEP), the Krieger-Li-Iafrate (KLI) and 

QO ' common energy-denominator approximations (CEDA) to the OEP, and the effective local 

potential (ELF). In particular, we show that the Slater, KLI, CEDA potentials and the ELF 
can all be defined as solutions to certain variational problems. We also provide a rigorous 
derivation of the integral OEP equation and establish the existence of a solution to a system 
of coupled nonlinear partial differential equations defining the Slater approximation to the 
Hartree-Fock equations. 
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1 Introduction 

The Hartree-Fock exchange operator associated with a first-order electron density matrix 7 G 
^ ; B^{W? X R3) is a Hilbert-Schmidt operator on defined by 

V0€i2(R3)^ (A»(r)=-/ 2im0(r')dr' (1.1) 

(see e.g. [3 [221 EH [32] for mathematically oriented introductions to the Hartree-Fock model). The 
construction of this nonlocal (integral) operator is the most computationally demanding step in 
solving the Hartree-Fock equations, especially for periodic systems. In the early days of quantum 
chemistry, Slater proposed to approximate the operator if by a local (i.e. multiplicative) operator 
Vxiy)., termed the exchange potential [31]. Nowadays, the nonlocality of the Hartree-Fock exchange 
operator is rarely seen as an obstacle for numerical calculations, at least in Gaussian basis sets. 
Recent advances in the development of orbital-dependent Kohn-Sham density functional [H [9l 
[TOl [131 [16] have caused a resurgence of interest in representing the exchange interaction by a local 
potential. Local exchange potentials have also been used as input in time-dependent and linear 
response methods for computing excitation energies and other properties [l4[ [T8[ [25} [35] . 

The first local approximation to the Hartree-Fock exchange operator originally proposed by 
Slater in 1951 [31] is given by 

^'.s(r)^-4-/ 4^dr', (1.2) 
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where p(r) = 7(r, r) is the electron density. The second approximation put forth in the same 
paper is 

^^.,xa(r) = -ap(r)i/3, (1.3) 

where is a positive constant. 

Motivated by Slater's work, Sharp and Horton [30] introduced a variational method for obtain- 
ing local potentials that approximate the Hartree-Fock exchange operator. Considering a local 
potential W and the associated one-electron Schrodinger operator Hw = — + W , they defined 
the energy functional E{W) as the Hartree-Fock energy of the Slater determinant constructed 
from the lowest eigenfunctions of Hw ■ The local potential W'^^^ which minimizes E{W) is 
called the optimized effective potential (OEP). The exchange part Wk.oep of W'^^^ obtained by 
subtracting from T/F*-*^^ the nuclear and electronic Coulomb potentials can then be used as an 
approximation to the Hartree-Fock exchange operator. This idea was elaborated by Talman and 
Shadwick [34]. Unfortunately, it is difficult to give a proper mathematical formulation of the 
(infinite-dimensional) OEP search problem or to build consistent finite-dimensional approxima- 
tions to the OEP (of course, the two issues are closely related). An alternative is to solve the 
Euler-Lagrange equation of the OEP minimization problem, the so-called integral OEP equation. 
Two well-known practical approximations to this equation are the Krieger-Li-Iafrate (KLI) ap- 
proximation [15] and the common energy-denominator approximation (CEDA) [8] which yield, 
respectively, the local exchange potentials Wk,kli and Waj.cEDA- 

More recently, several other techniques for generating local exchange potentials have been pro- 
posed im [33] . These techniques associate a density matrix 7 with an effective local potential 
(ELP) that provides a variational approximation to the corresponding nonlocal Hartree-Fock ex- 
change operator. The construction of the Slater potential or an ELP requires a density matrix 
7, whereas the OEP, KLI and CEDA potentials are obtained by solving a self-consistent system 
of equations in which the density matrix is unknown. It has been pointed out [HI [33] that the 
CEDA potential coincides with the potential obtained by iterating the ELP procedure until the 
orbitals and the ELP are consistent with each other. The CEDA potential and the self-consistent 
ELP also coincide with the localized Hartree-Fock (LHF) potential introduced in Ref. [29] 

The purpose of this article is to study the mathematical properties of the local exchange po- 
tentials described above. Let us emphasize that the mathematical problems under examination 
are nontrivial and that only partial results have been rigorously established so far. Although 
several articles containing mathematical statements about local potentials have been published, 
very few of those articles involve mathematically rigorous arguments, and some contain erroneous 
claims. Many current behefs about local exchange potentials rest on implicit assumptions (such as 
the existence and uniqueness of the solution, differentiability of a functional, existence of a Hmit, 
convergence of an asymptotic expansion, etc.) that cannot be or have not been formally proven. 
Our results may seem weaker than those presented elsewhere, but at least they are rigorous. 

The article is organized as follows. In Section [21 the main mathematical properties of the 
Hartree-Fock model are briefly reviewed. The Slater potential (|1.2p is dealt with in Section [31 
We provide a variational characterization of the Slater potential constructed from a given 7 and 
study its asymptotic behavior. Then we focus on self-consistent equations obtained by replacing 
the Hartree-Fock exchange operator with the Slater potential (|1.2p in the Hartree-Fock equations. 
These equations, first written by Slater [31], do not seem to be the Euler-Lagrange equations of 
any optimization problem. Therefore, we use a fixed-point method to prove that, in the radial 
case (one nucleus at the origin and spherically symmetric orbitals), they actually have a solution. 
Note that the situation is completely different if one uses potential l|1.3p instead of l|1.2p . In that 
case, the self-consistent equations have an unambiguous variational interpretation: they are the 
Kohn-Sham equations obtained with the exchange-correlation functional 

i?.,x.(p) = -^ / P^/^(r)dr, 

which has been extensively studied from a mathematical point of view in Ref. [3]. 
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In Section m we focus on the OEP. We summarize the known mathematical results for the OEP 
problem and provide a rigorous derivation of the OEP integral equation. We also study the KLI 
and CEDA potentials. In Section[5]we prove that the self-consistent ELP coincides with the CEDA 
potential. We do not provide a complete mathematical study of the self-consistent KLI, CEDA 
and ELP equations, but only examine the analytical properties of the corresponding exchange 
potentials. We prove that, given a set of molecular orbitals and under some technical assumptions 
(always satisfied in practice), the KLI potential and the ELP (hence the CEDA potential) are 
uniquely defined, up to an additive constant. 

In order to keep the notation as simple as possible, we focus in Sections [JMl on fully spin- 
polarized electron densities, i.e., systems with spin-up electrons only. With the notable exceptions 
of one-electron systems and two-electron triplet states, very few systems of practical interest are 
fully spin-polarized. However, the mathematical results stated in Sections [JUS] are completely 
general in the sense that they can be easily adapted to closed-shell and spin-polarized systems. 
Details are given in Section [6l 

All the proofs are postponed until Section [71 Basic concepts of functional analysis that are 
necessary to understand our arguments are summarized in the Appendix for the non-mathematical 
readership. 



2 Hartree-Fock exchange operator 

Let us first recall the mathematical formulation of the Hartree-Fock model. As we deal with fully 
spin-polarized systems here, the spin variable can be omitted. Without loss of generality, we will 
also assume real- valued functions. In the Hartree-Fock setting, the electronic state of a system of 
N electrons is described by a collection <I> = (?!>j)i<i<Ar of N L^-orthonormal orbitals: 

Mr)(f>^{r)dr = S,^, (2.1) 



or, equivalently, by the density matrix 

N 

7$(r,r') = 5]0.(r)^.(r'), 

1=1 

the electronic density being given by 

N 

p$(r)=7<,(r,r) = ^|0,(r)|2. (2.2) 

i=l 

Denoting by 

k—l 

the potential generated by the nuclei {zk is the charge of the fc-th nucleus, its position), the 
Hartree-Fock functional reads 

^"'(*) = IE/ \^Uv)?dv+[ Kuc(r)p*(r)dr+i / / P^^^}P±^ dv dv' 



2 ./b3 ./b3 r — r' 



Each term of the Hartree-Fock energy functional is well-defined provided $ G (if^(M'^))^, that is, 
provided V^i G {L'^{R^))^ for all 1 < i < TV, or in other words, provided the kinetic energy of $ 
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is finite. This property results from the inequalities [6] 

1/2 

I , (2.4) 



Jr3 |r-Hfc| xt^i-'^^ J 



and 

drdv'<l f P^^P^p-drdr' <N^^' ly f \WMr)\' dr] .(2.5) 
Jr^Jr^ \r-r'\ V^-'k^ / 



l7*(r,r')P 
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The Hartree-Fock ground state energy of the system is obtained by solving the minimization 
problem 

I^^ = inf {E^^{<S>), $ e Xn} (2.6) 

where 



The inequalities (|2.4p and l|2.5p imply that the Hartree-Fock functional is always bounded from 
below on Xpf. The Hartree-Fock ground state energy is therefore well-defined for any molecular 
system of arbitrary charge. The existence of a Hartree-Fock ground state, that is, of some <i>^® G 
satisfying 

has been proved by Lieb and Simon [22] for neutral systems and positive ions. It is also known 
that 'very negative' atomic ions are not stable: Denoting by Z the charge of the nucleus, 

• (12.61) has no minimizer when N > 2Z + 1 1211 : 



• there exists a universal constant Q > (whose optimal value is not known) such that for 
N > Z + Q, (1221) has no minimizer [32] . 

>From a physical point of view, the instability of very negative ions results from the fact that 
the excess electrons escape to infinity. Mathematically, it is due to a loss of compactness at 
infinity. These two viewpoints can be linked by the concentration-compactness theory due to P.L. 
Lions [23]. The existence of a Hartree-Fock ground state for 'moderately' negative ions is still an 
open problem (only the stability of H~ has been mathematically estabHshed) . 



The Euler-Lagrange equations associated with the minimization problem l|2.6p read 

N 

^<i>4>i = y^^ijlpj^ (2.7) 

J = l 

where A ~ i^ij) is a symmetric N x N matrix (it is the Lagrange multiplier of the matrix 
constraint (|2.ip ). and where J^<s, is the Fock operator 

1 1 

2 |r| 

In the above expression, ★ denotes the convolution product: 

(/*5)(r):= / /(r')g(r-r')dr', 

Jr3 

and K<i, is the so-called exact-exchange (or Hartree-Fock exchange) operator. It is the integral 
(nonlocal) operator defined by 



{K^<j,){v) = - f 2i(£lil)0(r')dr'. (2., 



4 



It is easy to see that K,s> is a self-adjoint Hilbert-Schmidt operator on L^(K^). Indeed, the kernel 
'''I'^^^.^l ^ is a square integrable function on x M.^ [22] . 

For neutral systems and positive ions, JF$ (for any <& e Ajv) is a self-adjoint operator on L^(M^) 
with domain and is bounded from below. Its essential spectrum is CTcss (•?^* ) = [0, +00). 

For positive ions, the discrete spectrum of J-c!> consists of an infinite non-decreasing sequence of 
negative eigenvalues of finite multiplicities, which converges to zero [24] , 

Any minimizer of (|2.6p satisfies the Euler-Lagrange equations (|2.7p . Using the invariance of 
the Hartree-Fock problem with respect to orbital rotation [27j, it is possible to diagonalize the 
matrix A = [Ay] appearing in (|2.7p . More precisely, if U is an orthogonal N x N matrix (i.e. 
such that U'^U = UU'^ = In) and if $ e Xn, then ${7 = iJ2^=lUj^(f>J)l<^<N e Xn and 
i?^^ (<!>[/) = £'^^($) (in fact 7^(7 = 7$, so that one also has T^u = •^$)- Let $ be a solution 
of i|2.7p and U an orthogonal N x N matrix which diagonaHzes the matrix A, i.e. such that 
U'^AU = Diag(77i, • • • , r]N). Then ^ = {ip,) i<i<N = is a critical point of (|2.7p . with the same 
energy as $, such that for all i, 

This means that 4" is a collection of TV orthonormal eigenvectors of the Fock operator. Besides, it 
can be proved that if $ is a Hartree-Fock ground state, then 

• Aufbau principle [22]: the rji& are the lowest N eigenvalues of (including multiplicities), 
i.e. rji = ti (up to renumbering of the orbitals); 

• No-unfilled-shell property [1]: cat < e^v+i, i.e. there is always a gap between the highest 
occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO). 

Consequently, solving the system 

^ = {(l}i)i<i<N e Xn, (2.9) 
ei < £2 < • ■ ■ < are the lowest eigenvalues of 

and applying orbital rotations to the so-obtained solutions provides all the global minimizers 
of l|2.6p . as well as, possibly, local minimizers and other kinds of critical points. 

Let us finally mention that it is possible to reformulate the Hartree-Fock problem is terms of 
density operators. Recall that the density operator T associated with the density matrix 7 is the 
self-adjoint operator defined by 

(T0)(r)= / 7(r,r')</'(r')dr'. 

In other words, the density matrix 7 is the kernel of the integral operator T. If $ e Xn, the 
density operator T$ associated with the density matrix 7$ is the rank-A^ orthogonal projector 

N 
i=l 

The Hartree-Fock energy functional can be written as a functional of the density operator: 

where 7 is the kernel of T and p^{v) — 7(r, r). If 7 is regular enough, 

Tr (-^At) = -1^^ Ar7(r,r')lr,.r ^r. 
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The above definition of Tr (— ^AT) can be extended to any non-negative self-adjoint operator T, 
by noting that 

Tr (-^AT^ = ^Tr (|V|T|V|) 

when 7 is regular, and since the operator |V|T|V| is self-adjoint and non-negative, the right-hand 
side can always be given a sense in M+ U {+00} (it equals the trace of |V|T|V| if this operator is 
trace-class, and takes the value -l-oo otherwise). 

The Hartree-Fock ground state energy and density matrices can be obtained by solving 

inf {£"''(T), T e Vn} (2.10) 

with 

■Pw = {T e 5(i2(R3)) I = T, Tr (-AT) < 00, Tr (T) = N} . 

A remarkable property of the Hartree-Fock functional [19] is that the minimizers of l|2.10p coincide 
with those of 

inf {£:"^(T), T e Vn} (2.11) 

where 

Vn = Convex hull ofVN = {T e S{L^iW^)) | < T < 1, Tr(-AT) < 00, Tr (T) = A^} . 
Recall that the notation < T < 1 means 

V0gL2(R3), 0< (</)|T|0) < ii^ii^,. 

Note that a generic element of Vn is of the form 

+00 

1=1 

where (?/>,) is a Hilbert basis of L'^{R^) with G ^^(M^)^ < n., < 1, and T,t^ = TV. 
This property is the theoretical foundation of efficient algorithms for solving the Hartree-Fock 
problem [4l[5]. 

In what follows, we will denote respectively by and the Hartree-Fock exchange operator 
and the Fock operator associated with the density matrix 7: 

(A'^0)(r) = - / 0(r') dr', = -^A + Kuc + P7 * A + 

Jr3 \r - r \ I |r| 



3 Slater exchange potential 

The Slater exchange potential associated with some <f> G Ajv has the following expression [31] : 

P*(r) iR3 |r-r'| ^ ' 

Obviously, the above definition does not make sense if p<i.(r) = 0. This is not a problem if > 
almost everywhere, since, in view of Proposition 13.11 below, l|3.ip defines an essentially bounded 
function on W" . If p$ vanishes on a set of positive measure, the Slater potential will be set 
to zero on fi. There is some arbitrariness here, but as the density of physical systems is positive 
almost everywhere, this is not an issue. 
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Note that in the case = 1 and > almost everywhere, the Slater potential is given by 



1 



\Mr)Mr')\ 



/\\2 



■dr' 



r - r' 



L(r' 



/M2 



3 r - r' 



dr', 



and therefore cancels out the Coulomb potential (this is a case of exact self-interaction correction) . 

The following Proposition collects the main mathematical properties of the Slater potential 
associated with a given <I> G Xn- 

Proposition 3.1. Let $ ~ {4>i)i<i<N G -^jv- 

(1) The Slater potential w* g is an essentially hounded function which satisfies almost everywhere 



Ur' <vls{r)<0. 



In particular, g vanishes at infinity. 



Besides, if p<i> > almost everywhere and if one of the conditions below is satisfied: 

• the orbitals 4>i are radial (i.e. spherically symmetric); 

• there exists l<p<3/2<q<2 .such that |r|p$ G LP{M.^) n 

the asymptotic behavior of the Slater potential is given by 



(3.2) 



(2) If pq, > almost everywhere, the Slater potential v'^ g is the unique minimizer of the variational 
problems 

inf {/|(w), w € i='(R3) +L~(IR3)} and ini {jI,{v), v £ L^{W^) + L°°im.^)} , (3.3) 
where 



Here and below, &2 denotes the vector space of the Hilbert-Schmidt operators on L^(M^) and 
Ilea the Hilbert-Schmidt norm (see Appendix). In particular 



. N r ,^ , I 7*(i",i"")7*(i"",r') 
w(r)7$(r,r ) + / ■ dr 



r - r" 



dr dr' , 



'(;(r)7$(r,r') + 



7<i>(r,r') 



drdr'. 



The condition that there exists l<p<3/2<q<2 such that |r|p$ G Lp{M.^) n L'>{ 
obviously holds true when p$ decays exponentially fast, which is the case, in particular, when $ is 
a Hartree-Fock |22j or a Kohn-Sham LDA ground state, or a solution to the self-consistent Slater 
equation p.4|) (this is a straightforward application of the maximum principle since Kohn-Sham 
LDA and Slater potentials vanish at infinity). 

In general, v'^ g is not a continuous function. This can be seen writing ^^^(r) as 



N 



<5(r)= - E 



r - r' 
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The functions 

|r — r'l 

are continuous, while the functions '^'p'"j('r)^''^ can be discontinuous at a given point, either be- 
cause one of the (jii is discontinuous, or because vanishes. It is worth mentioning two special 
cases. When $ = [(pi) is a Hartree-Fock ground state, it can be proved by elliptic regularity 
arguments (see [23 for instance) that for alH = 1, . . . , A^, 

(t>, e c°°(]r3 \ {R^,}) n c°'1(m3). 

It then follows that w*^ is globally Lipschitz in any compact set of 'M? \ /O$^(0) and C°° in 
M? \ (pJ^(O) U {Rfc}). Stronger regularity can be obtained if <f> = {4>i) is a Kohn-Sham LDA 
ground state, or a solution to the self-consistent Slater equation (|3.4p . In this case indeed p$ 
is positive in (the ground state of the corresponding mean-field Hamiltonian is positive and 
non-degenerate), so that v'^g is globally Lipschitz in and C°° away from the nuclei. 

We have not been able to recognize in the self-consistent Slater equations 
2 r 



0,(r)0,(r)dr = (5,„ (3-4) 
El < ■ ■ • < ejv are the lowest TV eigenvalues of (^\^ + ^mc + P* * + ^|,s ^ 



the Euler-Lagrange equations of some minimization problem. It is however possible to prove by a 
fixed point method that l|3.4p has at least one solution for neutral atoms and positively charged 
atomic ions, provided only radial orbitals are considered. Recall that a function <j) is said to 
be radial if there exists a function such that (/)(r) = <p(|r|). We will denote by L'^,{M?) (resp. 
Hl{m.^)) the set of radial L'^{m.^) (resp. radial H^{n^)) functions, and set 

Theorem 3.2. In the case of a single nucleus of charge Z > N , | fg.^[ ) has a solutio'^ <& = {(jii) G 
with EN > and the corresponding exchange potential ^ is globally Lipschitz in M.'^ , C°° 
away from the nucleus, and satisfies, for all ?7 > 0, 



<^(r) = -^+o(c-(^v-2^-'')M). 



Besides, the minimum of the Hartree-Fock energy over the set of the radial solutions to i3.4\ ) is 
attained. 

In summary, one can associate with any $ e Ajv a Slater potential v^g. Among all the local 
potentials that can be constructed in this way, we can select a few of them which might be more 
physically relevant than the others: 

(i) the potential v^^J , where <i>^^ is a Hartree-Fock ground state of the system; 

(ii) the potential v'^^g , where <i>^^ is a Kohn-Sham ground state of the system; 

(iii) the potential v^^g^ , where ^^'^^ is a solution to the self-consistent equations p.4p which 
minimizes the Hartree-Fock energy (over the set of the solutions to p.4p ). The existence of 
such ^^"-'^ is granted in the radial case for neutral atoms and positive atomic ions by the 
last assertion of Theorem 13.21 



^In the Aufbau condition "ei < ■ • ■ < ejv are the lowest N eigenvalues of (— ^ A + Vnuc + P* * + ""J g j"i the 
mean-field Hamiltonian is here considered as an operator on L^'' 
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4 Optimized Effective Potential (OEP) 



4.1 Original formulation of the OEP problem 

As already mentioned in the introduction, the OEP approach consists in minimizing the energy 
of the Slater determinant constructed with the lowest N eigenfunctions of some one-electron 
Schrodinger operator Hw = —^A + W, W being a "local potential". Note that in most articles 
deahng with OEP, the set of admissible local potentials is not defined. Leaving temporarily this 
issue aside, we denote by 3^ a given set of local potentials {y can be for instance the vector 
space y = L^^^{W^) + L°°(R'^) arising in the mathematical formulation of the density functional 
theory [20]). We introduce the set of admissible local potentials 



w= {w ey 



Hw ~ + is a self-adjoint operator on L^(]R^), bounded 

from below, with at least N eigenvalues below its essential spectrum j 

and the OEP minimization set 

X^^^ = 1$ = {0»}i<,<jv I (t>^ e H\M.^), (021) and (gSD hold for some W ew}, (4.1) 
where conditions 114.21) and 114.31) are defined as 

Hw4>i — ^i4'i: 

/ 0j(r)0j(r)dr = ^'^''^^ 
and 

ei < • • • < e^r are the lowest TV eigenvalues of Hw- (4-3) 
The optimized effective potential problem then reads 

/OEP = inf {i?«E($), $ e A-OEP} . (4.4) 

Denoting by $0^^ a minimizer to (|4.4p . an optimal effective potential is a local potential W'^^^ £ 
W which allows one to generate $0^^ through l(4!2l) - l(4!3| . 

In order to emphasize the mathematical issues arising from the above formulation of the OEP 
problem, it is worth recalling the general method for proving the existence of solutions to a 
minimization problem such as i|4.4p . The first step consists in considering a minimizing sequence, 
that is, a sequence (<i>")„gN of elements of X^^^ such that 

lim ^HF(^„) ^ inf ^ ^ ;^OEP| _ 

n — ' oc> 

It is easy to check that the sequence ((f>")„gN is bounded in (iJ-^(R'^))^, hence weakly converges, 
up to extraction, toward some <i>°° S {H^{M.^ j)^ . It is then standard to prove (see [22] for instance) 
that 

£;HF($oo) < jj^f {£;HF($), $ e X^^^} . (4.5) 

The difficult step of the proof is to show that £ X^^^ (if £ X^^^ , we can immediately 
conclude, using l|4.5p . that is a solution to l|4.4p ). For this purpose, we need to introduce a 
sequence (W„)„gN of admissible local potentials (W„ € W) such that $" can be generated by Wn 
via (|4.2p - (|4.3|) . If the set of local potentials y is e.g. a reflexive Banach space and if {Wn)neN 
is bounded in y, then {Wn)nGN converges (up to extraction and in some weak sense) to some 
potential Woo G 3^- We could then try to pass to the limit in the system 

Hw„€ = ere, 

C(r)0,"(r)dr = <5.„ 
e" < • • • < are the lowest N eigenvalues of Hw^ , 
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using more or less sophisticated functional analysis arguments, in order to prove that <f>°° satisfies 

— ' ' ' — ^Iv the lowest N eigenvalues of , 
hence belongs to X^^^ . 

To make this strategy of proof work, we therefore need to find a functional space y for which 
the sequences of local potentials generating the minimizing sequences of (|4.4p are bounded. Un- 
fortunately, we have not been able to find any non triviaH functional space y satisfying the above 
condition. This mathematical difficulty has numerical consequences: It is easy to construct dra- 
matic modifications of the (computed) optimized effective potential that are "almost solutions" of 
the OEP problem (make the potential oscillate and/or go to infinity at infinity), see e.g. [33] . 

4.2 A well-posed reformulation of the OEP problem 

A way to circumvent the issues raised in the previous Section is to replace (|4.2p - l|4.3p with formally 
equivalent conditions that do not explicitly refer to a local potential W [2]. 

Let us first deal with (|4.2p . For an operator W being considered as a local potential, it needs 
to be such that 

(W(l))^p = (W^)(j3 

for any regular functions (j) and t/i. This requirement is in fact, at least formally, a characterization 
of multiplication operators. Applied to Hw, this characterization reads 

iHw<l))4' - {Hwi'H = \ ('/'AV' - VA0) = idiv ((/)Vi^ - i^V0). (4.6) 

It is then clear that if $ = (00 e {H^{R^))^ satisfies (02]) with an operator Hw for which gH) 
holds true, we also have 

div (0iV0i - 0iV0i) = Ci(^i(/)i, 

(4 7) 

(j)i{r)(l)j{r) dr = (5^-, ^ ' ^ 

with Ci = 2{ei - ei). Conversely, if $ = {0^} G (iJ^(R^))^ satisfies l|4.7p . then, at least formally, 
<i> satisfies 114.21) with for instance 




N N 



W=— — , (4.8) 

2/9$ 

ei = 0, and = Ci/2 for 2 < i < TV. 

The idea then is to replace condition (|4.2p with the formally equivalent condition l|4.7p , which 
does not explicitly refer to any local potential. 

In order to account for condition l|4.3p in the same way, we remark that for any $ G Xn and 

all 1 < i < A^, 

e C^iR'), I f Mrf |VV'(r)|2 dr = {Hw - e,) 



^It is of course possible to construct finite-dimensional functional spaces y for which 1)4. 4p . with A*^^^ defined 
by l|4.ip . has a solution. Reducing artificially the class of admissible potentials is however not a very satisfactory 
way to tackle the OEP problem. 
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where C^(M'') is the set of compactly supported C°°(M'') functions. It follows from the above 
equality (see [2] for details) that conditions l|4.2p - l|4.3p are equivalent to 

j ^ 0i(r)(/)j(r) dr = 5^, 
vj e C^(R^), VI < i < TV - 1, 

+2(e,+i-ei) f / ^{vf<j>,{vfdv-Y^([ ^(r)0,(r)0, (r) dr 

Combining the above result with the formal equivalence between (|4.2p and l|4.7p with q = 2(6; — ei), 
it is natural to introduce the optimization problem 

/OEP ^ inf {/;HP($), CO e A^OEpj _ (4 9) 

where 



l<i<N 



, e //^(R^), / <i),(j)j = 3 = ci < C2 < • • • < CAT < oo, 



V2 < i < iV, div (0iV0i - 0iV0i) = Ci0i0,, VI < i < iV - 1, V?A e C^(K''^), 

2^ 



JR3 ~^ \JK3 / Wr3 \JRa 



We thus have eliminated any explicit reference to a local potential. Note that for any reasonable 
definition of y, it holds 

The connection between the original OEP problem l|4.4p and its reformulation (|4.9p can therefore 
be stated as follows: If <^^^^ = {(t)f^^)i<i<N is solution to (|4.9p . and if the reconstructed potential 

N N 
$:^?EP^^pEP^^^^|^pEP|2 

H/OEP ^ 2^1 (4_^^) 

2 PjOEP 

defines an element of W, then cOOEP jg solution to (|4.4p and M^o^p jg g^j^ optimized effective 
potential. 

It is proved in [2] that for a neutral or positively charged two-electron system, problem (|4.9p 
has at least one global minimizer (OCiep. Unfortunately, we have not been able to establish whether 
or not the reconstructed potential (|4.1ip is a well-defined function. 

Let us conclude this section by remarking that (|4.10p yields 

jHF ^ jOEP ^ jOEP 

A natural question is whether these inequalities are strict or large. We are only aware of two 
partial answers to this question: 

• it is proved in [2] that in the case of a single nucleus of charge Z >2 and TV = 2 electrons 
occupying radial orbitals, 

^HF < 70EP. 

• a formal perturbation argument is used in [12] to show that /^^ — /OEP j^j, non-interacting 
electrons and that /^^ <; jOEP -^^ ^j^g presence of an infinitesimal Coulomb repulsion term. 
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4.3 The OEP integral equation and its approximations 

The functional W ^ E^^{(j)f', • • • , 0]^), where (0f , • • • , 0]^) satisfy 

</.r(r)0f (r)dr = 5,„ (4.12) 

^ ■ ■ ■ 5: ejv the lowest N eigenvalues of H\y, 

is not well-defined for two reasons: First, the set of admissible local potentials has not been 
properly characterized, and second, l|4.12p may have several solutions if there is no gap between 
the HOMO and the LUMO of Hw- It is therefore a fortiori not possible to define the derivative 
of this functional. One can however give a rigorous meaning to the functional and its derivative 
for local potentials W satisfying the following assumption. 

Assumption 4.1. The potential W belongs to L'^{M.'^)+L°°(M.^), and the hamiltonian Hw ■ defined 
on the domain D{Hw) = -ff^(K'^), is a self-adjoint operator on L^(R'^), bounded from below, with 
at least N eigenvalues (including multiplicities) below its essential spectrum, and there is a gap 

V = - > (4.13) 

between (the N-th eigenvalue of Hw) and e]v+i (the {N + l)-th eigenvalue of Hw, or the 
bottom of the essential spectrum if Hw has only N eigenvalues below its essential spectrum). 

Under Assumption 14.11 the ground state density operator of Hw is uniquely defined and 
satisfies 

Tw - arginf{Tr(i/wT), TeT'Ar} 

= X{-oo,eF]{Hw), 

where X{-oo,eF] is the characteristic function of the range (— oo, e^^] and 

eF = ie^ + e^^,)/2. (4.14) 

In addition, if w e L^(R^)nL°°(M^) is such that < ri/2, then W + w still is a local potential 

for which Assumption 14. II holds true (this follows from the Courant-Fischer formula [26]). In this 
case, 

Tw+w = arginf{Tr(iJw+tuT), T e Pat} 

= X{-oo,eF]{Hw+w), 

with ep given by (|4.14p . It is therefore possible to define the functional 

w i-^ £^^{Tw+w) 

on the ball 

B,,/2 = {we L\R^) n L°°{R^), llwlUinL- < ??/2} . 
For W satisfying Assumption 1, one can also define the exchange part of the potential W as 

1 



r 



where jw is the kernel of Tw- It is easy to see that G L^(M'^) + L°°(M'^). We are now in 
position to state the main result of this section. 
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Theorem 4.2. Let W be a local potential such that Assu'mption \4- l\ holds true. Then, there exists 
a unique function eL^{R^)r] H^iR^) such that 

£:HF(TM/+.„)=fHF(Tw-)+ / ^?^(r)u;(r)dr + O(|lw|li,nL==0■ (4-15) 

In particular, the function w i-^ {T\y+w) is Frechet differentiable at w ~ Q. Denoting by 
BP{z) = (z — Hwr)^"^ the resolvent of H\y, by C a regular closed contour enclosing the lowest N 
eigenvalues of Hy^/ (see Figure 1), and by tiy(r,r') the kernel of the finite rank operator 

R\z){K^„ - v'^)R\z) dz, (4.16) 

it holds {r) = ti4/-(r,r). Let {4i^)i<i<N be a set of N orthonormal eigenvectors of Hw associ- 
ated with the lowest N eigenvalues eY < • ■ ■ < e]v of Hw- Then 

N 

e^(r) = 2 (r) [(1 - Tw)[eT - (1 - Tw)Hw{l - TwT^l - Th.)(X,„. - ^^f )0r](r). 

i=l 



c 











hhH — \ — h- 

i pW pW 
V ^1 ^N/ 


N+l 




Figure 1: Integration contour C in the complex plane 



Let us come back to the formulation l|4.4p of the OEP problem. If no artificial restriction on 
the set of admissible local potentials is enforced, y must be such that for all W ^ y, and all 
w € C^(]R^), one also has W -\- w £ y . Let us now consider a local potential W € y satisfying 
Assumption 14.11 Then W G W and it follows from (|4.15p that if W is an optimized effective 
potential, then 

= 0. (4.17) 

Although not obvious at first sight, l|4.17p is a rigorous formulation of the OEP integral equation 
introduced in |30} I34j. To clarify this point, we now assume that the spectrum of Hw is purely 
discrete (this assumption is impHcit made in [30l[34], but is obviously violated for isolated molecular 
systems, since for such systems, W is expected to have a constant, finite value at infinity). In this 
case, there exists a Hilbert basis )neN of L^(]R'^) consisting of eigenvectors of Hw associated 
with the eigenvalues < < • • • < e^ < e]v+i — ' ' ' ) the resolvent can be rewritten as 



Then, 



R'iz) = E 



[(1 - Tw)[e7 - (1 - 'rw)Hw{l - TwTHl Tw)(i^,„ - )</>r](r) 



^ I K^w - 

a=N+l » 



so that condition (|4.17p also reads 

i=l a=N+l 



0f (r)0r(r)=O. (4.18) 
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We have found no mathematical result on the existence and uniqueness of the solutions to the 
self-consistent OEP integral equation in the literature. 



The exact OEP integral equation is of little interest in itself, but can be useful to build ap- 
proximations of optimized effective potentials. It has indeed be proposed [SOj I34j . instead of 
approximating the optimization problem l|4.4p . to approximate the first order optimality condi- 
tion (|4.18p . This is the foundations of the KLI (Krieger-Li-Iafrate) and CEDA (common energy 
denominator) approximations. 

Starting from the idea of Sharp and Horton [30], Krieger, Li and lafrate proposed the following 
approximation of the OEP integral equation: 

N 

E E (0f iif*"--^'.*KLii0r)0r(r)0f(r)=o. (4.19) 

i=l j£N* , jT^i 

We refer to [15] for details on the derivation of equation (|4.19p . Note that the inner sum runs over 
all (occupied and virtual) states different from i, which implies in particular that equation l|4.19p 
is not invariant with respect to orbital rotations. In order to cover the general case when H\y has 
a non-empty continuous spectrum, it is preferable to rewrite equation l|4.19p as 

P*-(r)<KLi(r) = - / ^f^^rfr' + f:(0r|<KLi-i^*w|^r)l'^r(r)r (4.20) 

The KLI exchange potential w^kli can then be obtained by solving the self-consistent equations 

< • • • < e]v are the lowest N eigenvalues of H\y , 



P*w(r).fKLi(r) = - / '^T^'ji^'' + E<^^ ' <KLi - K^^^ I \<f>7ir)\ 



a.e. 



(4.21) 

Note that if ($'^ 

I ^a: KLi) solution to the above system, so is ($ j^'^kli ^) any real 
constant A. We have not been able to prove an existence result for (|4.2ip . Assuming that l|4.2ip 
has a solution ($^,f*KLi) ^ith <I>^ € Xn, it is however possible to prove the following: 

Proposition 4.3. Let ($^,w*kli) ^'^ « solution to I4.21\ ) such that G Xn and < 
min acssiHw)- Then p$ is a continuous, positive function on M"^, and i'*kli continuous, 
bounded function on M.'^ . Besides, 

(1) the potential i'^kli unique solution, up to an additive constant, to the minimization 
problem 

inf {J^l^\v), V e L\W^) + L°°{W^)] (4.22) 

where 

N 

In particular, $^ being given, the KLI potential is uniquely defined up to an additive constant; 

(2) it holds 

vtl^A^) = <I (r) + E (af - \K,^- I ^D) (4.23) 



"(«) = \{\\{v~ A^*)T*fe, - |(0,|(« - 



1=1 
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where u*^ is i/ie Slater potential associated with <& anrf where a = {af ) G R satisfies 

(/^-^*"')a*"'=/3*^^ (4.24) 

(3) the solutions of the linear system {4-24^ form a one- dimensional affine space of the form 

a* +M(1,--- ,1)^. 

Replacing a* with a* +A (1, • • • , 1)^ in \4-2S^ , amounts to replacing w*kli ^^'^ ''^'*kli + '^- 

Note that contrarily to the situation encountered with the Slater potential (see problem l|3.3p ). 
the quadratic functional J^^^ is convex (it is non-negative), but not strictly convex. A conse- 
quence of that is the non-uniqueness of v^^^, which is only defined up to an additive constant. 

Let us now turn to the CEDA potential introduced by Gritsenko and Baerends (3- This 
approximation consists in replacing in l|4.18p the denominators eY ~ ^Y with a constant, yielding 

N +00 

^ (</>ri^.»'-<cEDAka^>0r(r)0r(r)=o. 

i=l a=N+l 

Here also, it is possible to provide a more explicit formulation of this equation, still valid when 
Hw has a non-empty continuous spectrum: 

P*"- (l-)Wa;,CEDA(l-) = - / T 3] dv + {(j)^ \ W^^cEDA - K^^v \ (j)^ ) 0, (r) ■ (r). 

Jm \r-r\ 

(4.25) 

Let us also mention that the common denominator approximation amounts to replacing in l|4.16p 
the resolvent i?°(z) = (2 — Hw)^^ with the resolvent -Rceda(-^) = ^ H^p^^)^^ of the operator 

H^^^^ = eT^jir + e (1 — T,j,H ) , 

where e and e lay respectively inside and outside C. 

The potential v^^^^ solves the self-consistent equations 

0r(r)0f (r)dr = 4 



< • • • < e]v the lowest N eigenvalues of Hw , 

W = Kuc + P$w' * |7| + *^:c,CEDA: 

P*-(r).^*.CEDA(r) = - / '^7' dr- (^f | vfcEUA - K^^^ I ^f) '^r(r) (r) a.e. 

|r-r| ^^^.^^ 

To our knowledge, the question of existence and uniqueness of the solution to the above system is 
still open. 

It turns out that wJ.ceda coincides with the so-called local Hartree-Fock (LHF) exchange 
potential v^:^^ , obtained by Delia Salla and GorHng on the basis of completely different arguments 
(see [29] for details) . We will see in the next section that it also equals the self-consistent effective 
local potential f Jelp uM- 
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5 Effective Local Potential (ELP) 



The effective local potential associated with a given $ e was originally defined as the local 
potential minimizing the function |13) 



N +00 



1=1 a=N+l 

{4>a)a>N+i being a Hilbert basis of the orthogonal of the vector space generated by {4>i)i<i<N ■ A 
simple calculation shows that S'(i.(u) = J^^^{v) where 

[A, B] = AB — BA denoting the commutator of the operators A and B. An intrinsic formulation 
of the ELP problem therefore reads 

inf { J|^P(w), V e L^W") + (5.1) 
Proposition 5.1. Let $ = {(j)i)i<i<N G Xjs; . Any solution w*elp 1^5. 1\) satisfies 

P*(r)<ELp(r) = - / + E m<^i.M - {m^m Mr) </',(r) (5.2) 

and ifte symmetric matrix M* = [(</'i|'y*ELpl'?^i)] '■^ solution to the linear system 

(/-A*)A/* = G*, (5.3) 

f 0i(l")'/>j(l")</>fc(l")0i(r) , /■ $ / N, / 4* /i,|7^|,\ 

Besides, if the orbitals 4>i are continuous and if the open set M"^ \ pJ^(O) is connected, then the 
solutions to ^5.S^) form a one- dimensional affine set of the form 

so that w*ELP uniquely defined, up to an additive constant, on the set where p<s> > 0, and can 
be given arbitrary values on the set where p$ = 0. 

Comparing l|5.2p and l|4.25p , one immediately recognizes that the self-consistent CED A poten- 
tial and the self-consistent ELP potential defined by 

0r(r)</>f (r)dr = %, 

< ■ • • < e]v are the lowest N eigenvalues of Hw, 
1 

W = Kuc + P$w * 1^ + ""x.ELP^ 

P^wW^^.elpW^-/ ^1^'' dr' + I .f.E,p - K^^ I ) 0r(r) 0f (r) a.e. 

Jr3 |r-r| ^^^.^^ 

coincide. As already mentioned, we are not aware of a proof of existence of the solution to this 
system. We can however use Proposition 15.11 to show that if ^ has a solution ('i'^, t^*ELp) 

with e Xn and if ei < mmacss(Hw), then w*elp can be obtained from $^ by solving an 
optimization problem, which has a unique solution, up to an additive constant (the proof follows 
the same lines as the proof of Proposition 14.31 <j>Y then is a continuous, positive function on R^, 
which impHes that p$w is positive and that the above connectivity condition is satisfied). 
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6 Extensions to the General, Unrestricted and Restricted 
Hartree-Fock models 

In the general Hartree-Fock (GHF) model, each molecular spin-orbital (jii is a complex-valued 
functiorlfl with spin-up and spin-down components, i.e. (pi S L^(R'^,C^). The orthonormality 
constraint (j2.1|) is replaced with 



(^,(r) • (^j(r)* dv = 5., 



I] 1 



where 



The density-matrix can then be represented by a 2 x 2 hermitian matrix 



with 

N 

1=1 

and the electronic density p$ is the sum of its spin-up and spin-down components: 

= +pi(r), p^(r) = 7^'^ (r, r), p^(r) = 74^(r, r). 

The Hartree-Fock exchange operator associated with 7$ is the integral operator on L^(R'^,C^) 
defined by 

V</)eL2(R^C2), {K^cP){v)^~[ ^^^(r,r')-</>(r')dr', 

where • denotes the usual matrix- vector product, and local exchange potentials are 2x2 hermitian 
matrices of the form 

w-'-T(r) ti-^^(r) 

The variational definition of the Slater potential given in Section [3] provides a natural way to 
define a Slater potential for the GHF framework: It is the local potential v which minimizes the 
Hilbert-Schmidt norm of the operator {v — i^$)T<i,. A simple calculation leads to 

<s(r) = -i?*(r)-iS*(r)-^[S*(r)-i?*(r)-iS*(r)i?*(r)] (6.1) 

p$(r) 

= -S$(r)i?$(r)-i ^ [5$(r) - R^{r)E^{v)R^{v)-^] 



i;(r) 



where 



i?$(r) 74,(r,r) and S$(r) = / - — 5— -7$(r, r') • 7(r, r')* dr'. 

Jr3 |r-r'| 



Within the unrestricted Hartree-Fock (UHF) model, each molecular spin-orbital is (gener- 



. Denoting by A'^ (resp. Np) the number of spin-up (resp. spin-down) orbitals, and 



ally) chosen real-valued and either spin-up, i.e. (pi{v) = ( '^ij^^ \ or spin-down, i.e 0i(r) 







■'GHF models are of particular interest for systems subjected to magnetic fields; for such systems, complex-valued 
wavefunctions are needed. 
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ordering the spin-orbitals in such a way that the first Na of them are spin-up, the UHF density 
matrix reads 

7$ r,r = '* ^ ' ^, 

V 7i(r,r') 

with 

jlHry) = g0l(r),/-I(r), 7i'(r,r') = ^ W0]v„+, W- 

i=l z=l 

Likewise, the UHF exchange operator is diagonal: 
K, = ^^ll' ) with V0 G L'iR% (i^r'^)(r) = - |7---(r,jT ^^^,^ 

It is then easy to check that in the UHF setting, the general formula (|6.ip reduces to 

*.TT 



with 



-:rW = -T^/ rfr', p|(r)=7S^(r,r). 



One recovers in this way the spin-up and spin-down local potentials originally introduced by Slater 
in [31]. 

In closed-shell models, each molecular orbital G L^(R'^) is occupied by one spin-up and one 
spin-down electrons. Denoting by Np = N/2 the number of electron pairs, it holds 



7;T(r,r') = 7i'(r,r') = ^ 0,(r)0,(r), plir) = pi(r), v^^l^r) = ^'^-^W, 
1=1 

llHr,r') = 7i^(r,r') = 0, «*j-^(r) = v^^^^r) = 0. 

Proposition l3.lt Theorem I4.2[ Proposition I4.3[ and Proposition 15.11 apply mutatis mutandis to the 
RHF setting, as well as to the spin-up and spin-down components of the UHF exchange operator 
and local potentials. As outHned above from the Slater potential, the variational characterizations 
(|4.22p and l|5.ip of the KLI and ELP potentials can be used to defined KLI and ELP potentials 
in the GHF setting. 

7 Proofs of the main results 

Throughout this section, we denote by Br the open ball of R"^ of radius R centered at 0, i.e. 
Br= {r e M^ |r| < R} and by = \ Br. 

In order to simplify the notation, we adopt in this section the usual abuse of notation consisting 
in denoting by the same letter an integral operator and its kernel. 

7.1 Proof of Proposition [3TTt Properties of the Slater potential 

If follows from the Cauchy-Schwarz inequality that 

2 



f\\2 



l7*(r,r')| 



N 
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In the set where p$ > 0, one therefore has 



■dr' 



■dr' < 



p$(r) 7r3 |r-r'| ' " p^{r) 
In order to establish the decay property, we rewrite v^^ as 

N 



-^dr' = <5(r)<0. 



0,(r)0,(r) f 0.(rO0,(r') ^^, 



remark that 



0»(r)0j(r) 



r - r' 



< 1, 



and conclude using the following Lemma. 
Lemma 7.1. Let <i> = {4>i)i<i<N G 

V^J{V)= I 



r/ien Vij vanishes at infinity. Besides, if the 4>i are radial or if there exists l<p<3/2<g<2 
such that |r| e n L''{m?), then 



Proof of Lemma \7.1[ Let us denote by pij — (f>i4>j. By Sobolev embeddings, pij G L^(R'^)nL'^(M'^). 
For all i? > and all r g such that |r| > 2R, one has 



\V^,ir)\ < 



\^d.' 



|r'|<fl 



\r'\>R 



\^dr<l 
r' - R 



1 



It then follows from the Young inequality and the Lebesgue dominated convergence theorem that 



|PylXB|j 



1 



< 



XBi 



L' 
XBi 



\Pv\XBf,* — 
+ IIIPylXBJjIlLi 



XBf 



Therefore, Vij vanishes at infinity. 



The case of radial orbitals can be dealt with using the Gauss theorem, which provides the following 
expression for the potential Vij : 



Indeed, 



I 



R3 max(|r|, |r'|) |r| 



< rr / \p.Ar')\dr' 

|r| J|r'|>|r| 











^V~\ i P v + 


/ 




|r| J|r'|>|r| 


J\r'\>\r 



P-M^dv' 



(7.1) 



We conclude using Lebesgue dominated convergence theorem. 
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Let us now prove (|3.2p in the general case (non-radial orbitals), under the additional assumption 
that there exists 1 < p < 3/2 < 9 < 2 such that |r| \pij \ € LP(R3) n L'?(R3). For all r € R^ 



|r|T/„-(r) - 4 



|r| - |r - r' , , , 
-—p^j{v)dv 



3 r — r 



< 



It suffices to show that the right-hand side vanishes at infinity. For all i? > et all r G M'^ such 
that |r| > R{R+l), 



r - r' 



\v'\<R 



r — r' 



\v'\>R 



r - r' 



< 



1 

'r 



fi3XB% * ■ 



where /y (r) = |r| |pij(r)|. We then use the same argument as above: 



< 



< WUjXBf^WL^ 



-f Xb^ 

J^jXBf. * -J-J 
+ WfijXB^JlLP 



Li' 



XBf 



Lp' 



where p' = (1 - p^^y^ £ (3, +00] and q' = {1 - g"^)"^ G [2,3). The proof of Lemma O is 
complete. □ 

Let us now turn to the proof of the second assertion of Proposition 13.11 For all v S i'^(R'^) + 
L°°(R'^), the operator ^7$ is Hilbert-Schmidt. Indeed 

{vj^){r,r') =v{r)j^{r,r') G ^^(R^ x R^) 

since |7<i>(r,r')| < /9<i,(r) p<[.(r') with p<j> G ^^(R^) n L^(R.^). One can thus define on L^{M.^) + 



L°°(R3) the functional 
For all V and h in L3(k3) ^ ioo(-jg3)^ 



/ N / n 7*(r, I"') 
t.(r)7j.(r,r')+ ^ ' ,/ 



drdr'. 



J| (f + /i) = J| (i') + / f (r)p$ (r) + 



|7(r,r 



3 r — r' 



Therefore, all the local minima of are global, and they are characterized by the equation 



w(r)p$(r) + 



|7(r,r')|' 



3 r — r 



dv' = 0. 



If p$ > almost everywhere, the Slater potential is the unique solution to the above equation, 
and therefore the unique global minimizer of J^. 

The fact that the minimizers of and are the same comes from the fact that 7I — 7$ 
implies (-fir$7$, d7$)s2 = (-f'^*, '^7$)s2 ■ D 



7.2 Proof of Theorem 13.2b Self-consistent Slater equation 

The strategy of proof is based on a fixed-point argument. Notice that variational methods cannot 
be used since (|3.4p seems to have no variational interpretation. 
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For all ?7 > 0, we consider the problem 



(7-2) 



< are the lowest N eigenvalues of (-^A - ^ + p<i>v-ky^+ ""T^f^ (o" 

where 



' p$(r)+77 7R3 |r-r" 



The proof of existence of a solution to l|7.2p for ?/ = follows the lines of the proof of Theorem III. 3 
in [21]. We first construct, for 77 > 0, a continuous application whose fixed points are solutions 
to l|7.2p in X^. We then prove the existence of a fixed point of using Schauder Theorem. The 
existence of a solution to l|7.2p in the case when 77 = is finally obtained using some limiting 
procedure. Note that we have introduced the parameter 77 both in the nucleus-electron interaction 
and in the Slater potential. In the former term, 77 plays the same role as in [IJ (i.e. it enables us 
to control the decay of the orbitals at infinity) . The role of 77 in the latter term is to ensure the 
continuity of the nonlinear application T'' for 77 > 0. 



First step. Construction of the application T''. 
Let 7; > and 



if = |* = (V'.)i<«<A^e (i7,,i(M3)) 



N 



< I 



N 



In denoting the identity matrix of rank N. The semidefinite constraint [J^s i^ii'j] < In means 



Vx G M^, ( [ X^Xj < |x|^ 



It is easy to see that if is a nonempty, closed, bounded, convex subset of the Hilbert space 
{H^iM.^))^, containing X^^. For * e K, we denote by 7*(r,r') = = 
7*(r,r) and 

2 |r| |r| 

As the potential = + P* * j7[ + belongs to 

i'(R3)+i-(M3) ^ {ly I Ve > 0, 3{W2,Woo) € L''{R') x L°°(R3), ||W^oo||l- < e, W ^ W2 + W^} 

it is a compact perturbation of the kinetic energy operator. By Weyl Theorem [26], aessiF^) = 
o'css(— = [0,00). Besides, using Gauss theorem and the inequalities — -j^ < < v^'s < 

0, one has < V;^ < . Hence, 

G'^"^ ■■= -lA-^<n< 5'' -^A - ^. (7.3) 

As the hydrogen-Hke Hamiltonian , considered as an operator on L^(]R'^), has infinitely many 
negative eigenvalues, so does F^ (this is a straightforward consequence of Courant-Fischer min- 
max principle). Besides, the eigenvalues of the radial Schrodinger operator F^ being simple, the 
spectral problem 



ei < ■ ■ ■ < Eat are the lowest N eigenvalues of F^ (on L^(M'')), 
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has a unique solution $ = (0;) in XJ^ C -fsT up to the signs of the orbitals (j>i. We can therefore 
define a nonlinear application T'' from K to K which associates with any ^ G K the unique 
solution ^ = {(f)i) G XJ^ C K to (|7.2p . for which 0i > in a neighborhood of r = 0, for all 
I < i < N (by the strong maximum principle, 4>i cannot vanish on an open set of M.^). 



Second step. Existence of a solution to (|7.2p for 77 > 0. 

Using standard perturbation arguments, it is not difficult to prove that is continuous (for 
the norm topology). Let us prove that is compact. Let (^") be a bounded sequence 
in K, and let $" = T^'^". There is no restriction in assuming that {^") converges to some 
e {H^(M.^))^, weakly in (iJi(K^))^, strongly in (L^^iR^))^ and almost everywhere. This 
implies in particular that the sequence (p*n * ITy "I" bounded in L°° and converges almost 

everywhere to p\s/v * 171 + ^^''s^ when n goes to infinity. Using again (|7.3p and denoting by e" the 



«-th eigenvalue of , one obtains 
2 



. N ^ 1 /■ r 7+ 

2E(ll^'^"ll^^-2(^ + '?))'-2(Z + r;)2<^ / |VCP- / ^p^„ <^,»<0. 



Thus, for all 1 < i < the sequence ((/)")„gN. is uniformly bounded in H^{M.^) (independently 
of (^'")), and therefore converges, up to extraction, to some (ji/ G iJ^(M'^), weakly in H^{M.'^), 
strongly in ^^^^(R'^) and almost everywhere. Besides, using (|7.3p and Courant-Fischer formula, 
one obtains 

jZ + vf . „ . _rf_ 
- ' - 2^2- 

Up to extraction, (e") therefore converges to some G [~'^^T?^' ~^]- Next, by Kato inequal- 
ity [26], 

-A|CI < -sgn(0r)AC = 2(6r-V^*\)ICI 

< (7.4) 



As, moreover, (5*") and (<!>") are bounded for the norm topology, {V^„(f)") is bounded in 
L2(M3)^ so that (0f ) is bounded in H^{R^), hence in L°°{R^). Consequently, it follows from (fM)) 
and the maximum principle that there exists i5 > small enough and M > independent of i 
and n, such that 

\(f>^{r)\ < Me"('^"*')'""'. 

This implies that ((/)")„gN* converges (up to extraction) to (j)^ strongly in L2(]j3^^ particular, 
= £ '^A'- It is then possible to check, using the convergence of (^P") to 5''' and the 

convergence - up to extraction - of ($") to and of (e") to e^, that 

and next, using the positivity of p^n * jTy + s ^ ^^'^ Fatou lemma, that 
liminf-/ |VCI' = liminf2 / {V^^ - e^mf 



As on the other hand. 



/ IV.^I'P < liminf / |VCP, 
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(^'") converges to strongly in (77^ (R'^))^, which proves that T'' is compact. It then follows 
from Schauder fixed point theorem ^6] that has a fixed point G X"^, which is solution to 
(17:211. 



Third step. Existence of a solution to l|7.2p for 77 0. 

Let (77,1) be a sequence of positive real numbers converging to zero. As the sequence of correspond- 
ing fixed points ($''") is uniformly bounded in [H^iM?))^ and as — ^^X^-T^ < e^" < 0, there is 
no restriction in assuming that (<!)''") converges to some <I>* S [H^iM?)) , weakly in (_ff^(R^))^, 
strongly in (L[^q^(R^))^ and almost everywhere, and that (e^") converges to e* < 0. Besides, the 
sequence is bounded in {H'^{m.^))^ , hence in {L°°{R^))^ . 

Passing to the limit in the equation 0^" = e^" 0^" yields 



Vx,st 



Assume that /^s p$* < N . As 

^*-^-2^-H^^*''"*M' 

one has, using Courant-Fischer formula, and denoting by Xi{A) the i-th eigenvalue of A, 

t\ = Hm e^" 

n— !-4-oo 

= lim \i [F'l'i,^^ 

< hm A; - — + p$„„ ★ — 

\^ 2 |r| |r| 

1 A ^ 1 

-2^-M+^*'*W 



2 |r 

2z2 

It follows that for n large enough, the sequence (e^") is isolated from zero. As (<!)''") is bounded 
in (L°°(M^))^, we conclude, reasoning as above, that there exists M S M+ and a > such that 
for n large enough 

|0:'"(r)| < A/c-"l'-|. 

This implies that ($''") converges to $* € (iJ^(R^))^ strongly in (^^(M^))^, and consequently 
that /ju3 p$* = iV. We reach a contradiction. This means that Jj^a p$* = and therefore that 

This proves that are orthonormal eigenvectors of i^,^,. The fact that Ej < • ■ • < are the 
lowest eigenvalues of F^* follows from Courant-Fischer formula. 

In view of Proposition 13.11 the Slater potential 5 is equivalent to — at infinity. This proves 
that e* < • • • < < 0, from which it follows that the orbitals 0* enjoy exponential decay: For 
all ?7 > 0, there exists M G such that 

|</.^(r)| < A/e-(V^-''/3)l'-l. 

Using l|7.ip . one obtains 



<,(r) = -l + o(c-(^x/^-'')H 



Lastly, the same arguments as above can be used to prove that the minimum of the Hartree-Fock 
energy over the set of solutions to (|3.4p is attained. □ 
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7.3 Proof of Theorem 14.21 : OEP Integral equation 

Straightforward computations show that 

S^^i^fw+w) = S^^i^w) + Tr [Tj^i-yw+w - Iw)) + ai-^w+w - "tw,"tw+w - 7vy), (7.5) 
where 

a{j„j,)^'-f [ 7i(r,r)72(r-,rO-7i(r,r072(r,r-) 
2./n3./B3 r-r' 



The second term of the right-side of (|7.5p is well-defined since both jw+w and are finite rank 

-1 



operators with range in i/^(R'^) = D{T-y^,) ~ Let 



6 = ( max||i?°(z) 

Denoting by = one has for all u) e B,^/2 = {w £ L^{R^) n L°°{m.^), HwlUinL- < ^/2} 

such that ||it;||i=o < S, 

R^{z) ^ {z~Hw+w)~^ = {z-Hw~w)-^ ^ {{z-Hw){l-R°{z)w))-^ = {1 - R°{z)w)-^R°{z) 
and 

(1 - R°{z)w)-^ - 1 = (1 - R°{z)wy^R°{z)w = R°{z)w{l - R°{z)w)-\ 
Using the complex-plane integral representation 

iw+w = R'"{z)dz, 



one is led to 

IW+w — IW 

Hence, 



= — (({R'"{z)-R'^{z))dz^— {{1- R'^{z)wy^R'^{z)wR'^{z)dz 
2m Jc 27ri Jc 

= — I R°{z)wR°{z) dz + — (f R°{z)w{l - R"{z)wy^R°{z)wR"{z) dz. 
27ri Jc 2tt\ Jc 



•^7w [iw+w -iw) = ^ ^ T^^^R^{z)wR''[z) dz + ^^£ :F^„R\z)w{1 ~ R\z)w)-^ R''{z)wR\z) dz 

^ -^T (f HwR^{z)wR"{z) dz + ^ {{K-rw - v^)R^{z)wR°{z)dz 
27ri Jc 27ri Jc 



ic 

+— / T-y„R°{z)w(l - R°(z)wy^R°(z)wR°(z) dz 
2tt\ Jc 



(f {-1 + zR°{z))wR^\z) dz + ^ (fiK^„ - vY)R°{z)wR°{z)dz 
27ri Jc 27ri Jc 

-h-^-T / J'-y„R°{z)w{l - R°{z)w)-^R°{z)wR°{z) dz 



= -WJW + — i zR"(z)wR°(z)dz + — (f(K^„ -vY)R°(z)wR°(z)dz 
27ri Jc 27ri Jc 

+ — I J^^,,.R°{z)w{l - R°{z)w)-^R°{z)wR°{z) dz. 
27ri Jc 

To proceed further, we make use of the following technical Lemmas, whose proofs are postponed 
until the end of the present section. 
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Lemma 7.2. For all z e p{Hw) , (1-A)i?°(z) and - A) are bounded operators on L"^ (M.^) 

and (1 — IS)EP{z) is the adjoint of R^{z){\ — A). Besides the functions 

z^ R"{z){l- A) and z (1 - A)i?°(z) 

are analytic from p{Hw) into C{Lr'{M?)). 

Lemma 7.3. 

(1) For all V G L^(M^), the operator (1 — A)~-'^u(l — A)"'^ is trace-class and 

||(l-A)-l^;(l-A^l|le, <^||«|Ui. 

(2) For all v € L'^{M.^), the operator v{l — A)~^ and its adjoint (1 — A)~^v are Hilbert- Schmidt 
and 

||«(l-A)-i||e.==||(l-A)-i.||e.- """"^^ 



(87r)i/2- 

Using the above two Lemmas, it follows 

Tr (j) zR°{z)wR"{z)dz^ = ^ (j) zTt {R°{z)wR°{z)) dz 

= / zTr (i?°(z)(l - A)(l - A)-^w{l - A)-\l - A)R°{z)) dz 

= ^ (f zTt ((1 - A)i?°(z)2(l - A)(l - A)-^w{l - A)"i) dz 
= Tr (^{1 - A) / zR°{z)^ dz) (1 - A)(l - A)-^w{l - A)-^ 

Denoting by 



2711 JC 



Hw^ / XdPx 
the spectral decomposition of ffvK, it holds 

1 /■ 1 /■ / /■+°° z 



■ dz dPx 



2t:\ Jc {z-\Y 
dP\ = 7w- 



Hence, 



(2^ jc '^^ = " ^)7w(l - A)(l - A)-^w{l - A)-i) = Tr {^iww). 



Tr 

We thus obtain 



+Tr ^ J^^„-i?°(2)w(l - R"{z)w)-^R°{z)wR°{z) dz^ (7.6) 

+Q;(7iv+u) — 7VK, 7VK+U) — 7vy)- 
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Let us denote by 



/3 = max||(l- A)i?°(z)|| 



AsWe L^{R^) + L°°{R^), one also has e 
L°°(M3) such that = V2+Voo. Then, 



Tr ( (K^^-v^)— j> R%z)wR%z)dz 



< 



i°°(R3). Let V2 e L2(R3) and £ 



Tr iK^„ - v^)— i R%z)wR^iz)dz 



27ri 



1 



with 



Tr (^{K^,,. - t;^)^ j( R°{z)wR\z) dz 



Tr ( V2^^j R^{z)wR^{z)dz 



< B^{\\K,„\\ + |koo||L~)ll(l - A)-V,(l - A)-i|le, 



27r 

< C\\w\\li <C\\w\\LinL^ 



(7.7) 



and 



Tr (^V2^^£R°{z)wR°{z)dz^ = Tr (^^^ j( ^2(1 - A)-i(l - A)ii'°(z)u;(l - A)-i(l - A)ii'"(z) dz^ 



< 



2tt 



\\v2{l-A)-'\\eAw{l~Ay'\\e 
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< C\\w\\l2 <C\\w\\LinL^ 



w ^ Tr ( {K^,,, - Vx)^^ f R'^{z)wR'^{z) dz 



(7.8) 



The Hnear form 



therefore is continuous on L^(]R'^) n L°°{M?). It remains to prove that the last two terms of the 
right-hand side of l|7.6p are 0(||ii;j||ip^^). The first one is easy to deal with. Indeed, 



Tr 



j> ^iwR^{z)w{l - R\z)wy^R^{z)wR\z)d:^ 



< 



|C|/33|j^,„,.(l_A)-i|| 



27r 1 



\w\\l^ |l(l-A)-Ml-A)-^||e, < 



|C|/33||.F,,„,(l-A)-i| 



167r2 ( 1 



The second term can be split as 

where D{-, •) denotes, as usual, the Coulomb energy 

/(r).g(r') 



\hw+io -jw){r,r')\ 



'\\2 



■ drdr', 



r — r' 



Dif,g) 



dr dr' , 



for which [26] 



3C e M+ s.t. V/ e L'/^R'), < D(/, /) < C\\f\\ls/r,. 
As both p-y„_^_^ and belong to L1(K3) n L'^{M.^), 

' PlW ' PlW + m PlW ) < c\\ 

^ ^\\P'yw+^ ^ PiwWl^ \\P'rw+^ ~ PlwWh'i ■ 
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We now make use of the following characterization of the norm [17], which is valid for all 

I < p < +00: 



geLp'i 



sup 



fg, 



where - + -^7 = 1. In our case, one obtains 



= sup 

96i~, ll<7lll-°°=l 

= sup 
< sup 

96i=", ilslll.°°=l 

= sup 

9ei'=°, llsllt°°=i 



Tr {{iw+w - lw)g) 
\\{lw+w - 7iv)g||ei 



/(I - i?"(z)u;)-ii?0(z)wi?"(z)5dz 
27ri Jc 



Si 



< 



167r2 ( 



1 - 



and 



l|P7 



sup / {Piw+^- Piw)g 

96-£'Ml9lll,2=l. ^ 

= sup Tr ((7vv+„ - 7H/)5) 

< sup \\{iw+w~iw)g\\ei 

96i^ Il9lli,2=l 

1 



sup 

ffei-^ Il9llr2=l 



sup 



27ri jc 



(1 - R°{z)w)~^R°{z)wR°{z)gdz 



Si 



< 



ICI/32 



||(l-A)-i,«i|e. ||(l-A)-ig|| 
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167r2 [l 



s 



Hence, 



< D{p^„^^ - P-twyP-rw+u. ~ Piw) < 



C 



Lastly, one obtains, using again Cauchy-Schwarz and Hardy inequalities. 



As 



liTw+io - 71^1162 < 1171^+10 - 71^1101 < 



167r2 1 



27 



and 



V-^ i R"{z)wR"{z)dz 
27ri J(, 



27ri 



J-^ I V(l - A)-i(l - l^)R°{z)w{l - A)-i(l - A)i?"(z) 
27ri 

— / V(l - A)-i(l - A)ii'"(z)w(l - i?"(z)u;)-ii?0(z)w(l - A)-i(l - A)i?°(z) dz 



< C i|u;||.2 + 



we conclude that 



\{lw+w -7w)(r,r')l' 



dvdv' = 0{\\w\\l^^^. 



We have therefore estabhshed that the Frechet derivative of the function w i— > £^'^ {'^w+w) is the 
Unear form 



R°{z)wR°{z)dz 



It follows from ([L7|)-(ILH1) that this Hnear form is in fact continuous on Li(M^)nL2(M^). Therefore, 
there exists e {L^{m.^)f^L'^{R^))' ^ L2(R3) + L°°(R3), such that for all w e L^{M?)r^L°° {M!^), 



Tr [{K,„ -v^)l^J^R\z)wR\z)dz^ = 



Q W. 



Using [i?°(z), 7via] = 0, the analyticity of the function z i~> (1 — jw)R'^iz){l — jw) in the interior 
domain defined by C, and Cauchy's formula [28], it is easy to show that 



1 
27ri 



R"{z)wR^{z)dz = 



27ri 



jwR°{zhww{l - 7H')i?°(^)(l - iw) dz 



27ri 



The left-hand side of the above equation therefore defines a finite-rank operator. Let {(pi^)i<i<N be 
a set of N orthonormal eigenvectors of Hw associated with the lowest N eigenvalues e < • • • < eJv 
of Hw ■ It holds 

N 

l,lR\z)wR\z)dz = Y,{\4,Y){cj^Y\)^^ i ^^{l-lw)R\z){l-^w)dz 
27ri Jc j-^ 27ri Jc z ~ ef 



1=1 

AT 



E(- 

^ \ 27ri 



c z-e 



w 



(1 - iw)R''{z){l - 7H^) dz I) 



Using again the analyticity of the function z i-^ (1 — ^w)R^{z){l — ^w) in the interior domain 
defined by C, and Cauchy's formula, we then obtain 



27ri 



1 



c 



(1 - 7vv)i?"(z)(l - 7vv) dz={l- -iw)[eY - (1 - 7H')-ffw(l - IwT^l - Iw) 



Multiplying the above equality by {K^„ — v^) on the left-hand side and taking the trace, we are 
led to 

Tr [^K,,,. - «f )^ j( R\z)wR\z) dz 



w 
g w 
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with 

N 



g'^ir) ^ 2 J2 '/'fW [(1 - lw)[eV - (1 - 7w)i?w(l - 7w)]~^(l - lw)(.K,,, ~ v^)cj^T]{r). 
i=i 

As the s are in H'^{R'^) and as the range of the operator [e]^ — (1 ~ lw)Hw{^ — 7iv)]~^ is 
contained in i?^(M'^), the function pw belongs to L^(M^) n i7^(M^). Using similar arguments, one 
can easily show that the operator Tw defined by l|4.16p is finite-rank and that pw(j) = twij^ii^)- 

It remains to prove Lemmas 17.21 and 17.31 

Proof of Lemma\2M Let z be in the resolvent set p{Hw) of Hw- By Assumption 14.11 D{Hw) = 
H'^{M?). Hence, {z — Hw), considered as an operator from 77^(]R^) to L^(]R^), is invertible. As 
W G L^(M'^) + L°°(R^), it is also continuous, hence bicontinuous in view of the inverse mapping 
theorem [26]. As so is (1 — A), (1 — A)_R°(z) is a bounded operator on L^(M^). 

On the other hand, it holds, for all c > such that (z — c) € p{Hw), 

i?0(2)(l - A) = R°{z){{z - c) - Hw)R\z " c)(c - A/2)(c - A/2)-i(l - A). 

The operators BP{z){{z — c) — Hw) = 1 — cBP{z) and (c — A/2)^^(l — A) are bounded operators 
on L2(k3)_ Besides, 

(c - A/2)-i(z -c-Hw) = -{l-{c~ A/2)-\z - W)) . 

As e L2(M3) + L°°(R3), one can write W as W = W2 + Woo with W2 G L^{R^) and Woo € 
L°°{E.^). The operator {c — A/2)~^{z — Woo) is a bounded operator and its norm goes to zero when 
c goes to +00. Lastly, the operator (c — A/2)~^W2 is Hilbert-Schmidt, and its Hilbert-Schmidt 
norm 

||(c-A/2)-H^.||,^.-^^^^drj mU., 

hence its norm in £(L^(M'^)), go to zero when c goes to infinity. The operator (c — A/2)~^(z — c — 
Hw) is therefore bounded on L^(M^) and invertible for c large enough. Its inverse, i?°(z — c)(c — 
A/2) also defines a bounded operator. This proves that i?°(z)(l — A) is a bounded operator. 

The analyticity of the functions z i-^ (1 — A)i?"(z) and z 1-^ _R°(z)(l — A) follows from the 

analyticity of the resolvent on the resolvent set: For zq G p{Hw) and z e p{Hw) such that 
|z-zo| < \\R°izo)\\-\ it holds 

+ 00 



i?°(z) = ^(z-zo)"i?°(zo)"+i. 

71=0 

□ 

Proo/ of Lemma Vr^ Let us first prove the second assertion. The kernel of the operator u(l — A)~^ 
is explicit and reads 

e-|r-r'| 



fc(r, r') = v{t)- 



47r|r - 

As 



- A)-i is Hilbert-Schmidt and ||u(l - A)-'^\\e.^ 



;87r)i/2 ■ 

In order to prove the second assertion, we write w as u = w+ — w_ with = max(w,0) and 
V- = max(— WjO), and introduce the operators A± — y/v±{l — A)^^. As y/v± £ L^{R^), the 
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operators A± are Hilbert-Schmidt and such that ||A±||e2 
A)-iu(l - A)"i = ^+^+ - A*_A^ is trace-class and 



\il-A)-\{l-A)-'\\< 



8tt 



V+\\l2 



-III 



t}I|U./(87r)i/2. Hence, (1 



8tt 



□ 



7.4 Proof of Proposition 14.3b Properties of the KLI potential 

As D{Hw) = H'^{R^), the eigenfunctions 4>Y are in i/^(R^), and are therefore continuous on R'^. 
Under the assumption that ej^ < min(Jcss{Hw), the ground state (/)]|^ is non-degenerate, and 
positive on M? . Consequently, p$ is continuous, and positive on M? , so that 



<KLi(r) = vf^ (r) + {i€ Kkli I €) ~ i^Y \K^« I €)) (7-9) 



is a continuous, bounded function on 



i=l 

3 



Proceeding as in Section mi one can show that the functional J^^^ is well-defined on L'^(IR^) -|- 
L°°iM?) and that the global minimizers v of l|4.22p are exactly the solutions to the KLI equation 



p^w{v)v{v) 



I 'T^'Ji^'' d-' + E(^r I V ~ I i^rwr (7.10) 



It remains to prove that the set of solutions to the above equation is a one-dimensional affine 
space. To this end, we remark that a potential 



N 



p$w(r) 



is solution to (|7.10p if and only is the vector a*"' = (af"') G is solution to the linear 
system (0211). We therefore have to show that Ker(/Ar - 5*''') = M(l, • • • ,1)'^ and that S 
Ran(/7v-S'*). 

Let 2/ e M^. One has 



N \2 / N \ / N 



■-ii.-s--)y = j:y^- H -^T.y^- 



2/ _ 

i=l i=l 

with equality if and only if, for all r G R^, there exists A(r) such that yicj)^ {r) ~ A(r)(^J^(r) for 
all 1 < i < A^. As (jj^ > on R^, this condition is equivalent to y = (j/^) S M(l, • • • , 1)^. Thus, 
Ker(7Ar — S ) = R(l, • • • , 1) . Lastly, using the fact that S is symmetric, one obtains 

Ran(/Ar - S"*"') = Ker(/7v - S'^"' )^ = |z = (z^) e M^, ^ 2, = o| . 
It is easy to check that /J*"^ e Ran(/Ar — S"**^). 
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7.5 Proof of Proposition 15.1^ Properties of the ELP 

For all V G L^{R^) + L°°{M.^), the operator B'^v = [v,-/^] is Hilbert-Schmidt. One can therefore 
define on L^{M.^) + L°°{m.^) the functional 

- ^ll[«-i^*,7*]ll|, = [A*,7*]ll|,- 
For all V and h in L^{R^) + L°°{R^), 



and 



p<i,(r)u(r) + 



'2 



The global minimizers v of l|5.ip are therefore exactly the solutions to the equation 



l7*(r,r')|^„, 



N 



(7.11) 



Multiplying the above equation by and integrating over M'^, one then observes that a function 
V satisfying 

p^{v)v{v) =~ [ dr' + f2 (M., ~ (0.|i^*|0,))0.(r)0,(r) 



is solution to l|7.1ip if and only if the matrix M is solution to the linear system 

{I-A^)M = G*. 



(7.12) 



Let us now prove that, if the orbitals are continuous and if R^ \ p^^(O) is connected, then 
Ker(/-A*) = R/w and G* € Ran(/-A*). For this purpose, let us consider a matrix M e Ms{T^) 
such that (/ — A*) M — 0. As M is symmetric, it can be diagonalized in an orthonormal basis 
set as 

M = Diag(Ai,- • ■ ,\n)V 

where J7 is a unitary matrix. Denoting by (V'l, ■ • ■ , V'Jv)^ = . . . , (jiw)'^ , a simple calculation 

leads to 



TV 

0=((/-A*)M,M)f = ^A?- / 



N 



dr 



where {■,-)f is the Frobenius inner product on Ms{N). As J7 is a unitary transform, the ipi 

N 

are orthonormal for the L^(M'^) inner product and ^^-04(1")^ = P*(i")- Therefore, using Cauchy- 
Schwarz inequality, 

N 2 



J2Mr?] EAfV..(r)M =p*(r)EA?V',(r)^ 

i=l I \i=\ j 1=1 
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with equality if and only if there exists C(r) such that Xitpi{r) = C{r)tpi{r) for all 1 < i < 
Hence, 



N 



N 



2 



N „ N 



with equality if and only if for almost all r e M"^, there exists C(r) such that Xiipi{r) — C{r)^i{r) 
for all 1 < i < TV. 

If the orbitals (jii are continuous, so are the functions ^pi. Let us consider the open sets fli = 
\ ij}-\0) and n = Uilif^^ = \ p-i(O). On Qi, one has C(r) = X,. This implies that the 
function C(r) is constant on each connected component of fJ. If ft is connected, one therefore has 
Ai = A2 = • • • = Atv, i.e. Af is proportional to the identity matrix. 

In summary, under the assumptions that the orbitals are continuous and that \ pJ^(O) 
is connected, 

(1) the linear equation (|7.12p has a solution if and only if G* S Ran (/ — A*). Note that 
Ran (J - A*) = Ker (j - (^*)*)"^ = Ker (/ - A*)^, since ^* is self-adjoint for the Probe- 
nius inner product. It then follows that Ran (/ — A*) = Span(/jv)^. Since {In,G'^)f = 
Tr(G*) = 0, G* e Ran (/ - A*) and 1(7121) has at least one solution Aff ; 

(2) if Mf is a solution to ([7A21) . then the set of the solutions of l(7?T2l) is {Mf + A/r« , A G R} . 
Note that replacing Af* with Af* + A/rn in (|7.12p amounts to replacing w*elp '^ith w*elp + 



Appendix: Elements of functional analysis 

This Appendix summarizes the basic concepts of functional analysis needed to understand the 
statements of the results contained in the present article. The additional concepts and results 
used in the proofs can be found in Ref. |26j . 

Let us first recall the definition of the functional spaces used throughout this article. In the 
following, all the considered functions are real- valued Lebesgue measurable functions on M"^. As 
usual, two functions which differ only on a set of measure zero are identified. 

For 1 < p < 00, the space is defined as 



Lf(R3) = (u / \u{r)\Pdr<oo 



Endowed with the norm 

L^iM?) is a Banach space. The space L'^iM?) is a Hilbert space for the inner product 

{u\v) = / u{y) v{y) dr. 



The space L°°(R'^) is the vector space of essentially bounded functions. A measurable function u 
is essentially bounded if there exists a constant M such that |u| < Af almost everywhere (a.e.), 
i.e. everywhere except, possibly, on a set of measure zero. Endowed with the norm 

|1u|1lo= =inf {A/ > I |u| < A/ a.e.}, 

L°°(R^) is a Banach space. One has for all u e L°°(M^), 

< a.e. 
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Forall 1 < p < g < oo, the space LP(R3)nL'?(M3), endowed with the norm H-liiPni^ = I1-|1lp + |M1l<j, 
is a Banach space. Likewise, for all 1 < p < g < oo the space 

LP{R^) + L^R^) = {u I 3(up, Uq) e LP(R3) x L'^(R^), u = Up + Uq} , 

endowed with the norm 

\\u\\lp+li = inf {llMplliP + {up,Ug) e LP{M.^) x L\S.^), u = Up + Uq] 

is a Banach space. 

In quantum mechanics, the kinetic energy of a one-particle wavefunction (j) is \ j^^ |V(/)p. It 
is therefore natural to introduce the vector space 

H\W^) = {u e L^{W^) I Vu G (i^(]R^))^} . 

Endowed with the inner product 



/ u(r) v{r) dr + Vw(r) • Vw(r) dr, 



H^{M.^) is a Hilbert space. We will also use the Hilbert space 



whose inner product is 



{u,v)h^ ^ I w(r)u(r)(ir+ / Vti(r) • Vi;(r) dr + / T i^) ^ (r) dr. 

The functional spaces i7^(R'^) and H'^{M?) belong to the class of Sobolev spaces. 

Lastly, LfQj,(M'^) is the vector space of the functions u such that |u(r)|^' dr < oo for all 
compact sets K cM."^. 

The second part of this Appendix is devoted to linear operators on L^(R'^). The set of the 
continuous linear operators from L^(R^) to L^(R'^), also called bounded operators on L'^{R^), is 
denoted by £(L^(M^)). The adjoint of a continuous linear operator T e is the unique 

operator of £(L^(]R^)), denoted by T* , defined by 

V(w, v) e L2(R3) X L^{W^), {T*u\v) = {u\Tv). 

The operator T € £(L^(R^)) is called self-adjoint if T* = T. The vector space of self-adjoint 
continuous Hnear operators on L^(R^) is denoted by 5(L^(R'^)). If T is a self-adjoint operator, it 
is usual to denote by 

{u\T\v) ^ {Tu\v) = {u\Tv). 
Let T G C{L'^{R^)) and (e„)„eN a Hilbert basis of L'^{M.^). The value of the sum 

is independent of the choice of the Hilbert basis (e„)„gN. The operator T is called Hilbert-Schmidt 

....... r....,S.. 



VnSN / 
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The set of Hilbert-Schmidt operators forms a vector space, denoted by 62. It is in fact a Hilbert 
space for the inner product 

{S,T)e2 ^(S'e„|Te„). 

The norm associated with (•,-)s2 is denoted by || • ||e2- It can be proved that T e C{L^{R^)) 
is Hilbert-Schmidt if and only if there exists a function of L^(R'^ x M.^), called the kernel of the 
operator T and usually denoted by T as well, such that 

{Tu){r) = f T{v,v')u{v')dr'. 

It holds 

\\T\\e,= (f \T{r,r')\^drdr' 

and T is self-adjoint if and only if T{r' , r) = T(r, r'). 

Let now T G S{L'^{M.'^)) be a non-negative self-adjoint operator (i.e. {u\T\u) > for all 
u e i^(M^)) and (e„)„gN a Hilbert basis of L'^{R^). The value of the sum 

^{en\T\en) 

neN 

does not depend on the choice of the Hilbert basis (e„)„gN- If this sum is finite, T is called 
trace-class and the trace of T is defined as 

Tr(r) = ^(e„|T|e„). 

riGN 

A (non-necessarily self-adjoint) operator T G £(i^(R'^)) is called trace-class if the non-negative 
self-adjoint continuous operator |r| ~ (T*Ty^^ is trace-class (the square root of a non-negative 
self-adjoint operator is defined below). The set of trace-class operators on L^(R'^) forms a vector 
subspace of 62, denoted by 61. Endowed with the norm 

||T||e, =Tr(|r|), 

61 is a Banach space. For all T e Si, the sum J2neNi^n\Ten) is finite and independent of the 
choice of the Hilbert basis (e„)„eN. The trace Tr defines a continuous linear form on Si. 

Most linear operators arising in quantum mechanics are not continuous linear operators. An 
example is the one-particle kinetic energy operator Tk = — |A. As the Laplacien of a function of 
L^(]R'^) is not, in general, a function of L^(]R^), Tk cannot be defined as a linear application from 
L^(]R'^) to itself. The useful definition of linear operators is the following: A linear operator T on 
L2(M3) is a L2(K3)..^,g^iyg^ ijjjgj^j. application defined on a subspace D(T) of L'^{R^). The set D{T) 
is called the domain of the linear operator T. For instance Tk is a linear operator on L^(R^) with 
domain D{Tk) = H'^{R^) (for all u e H'^R^), Au € L'^{R^), and Tku therefore is a function of 

L2(M3)). 

Let T be a linear operator on L^(R'^) with dense domain D(T). The adjoint of T is the unique 
linear operator on L'^iR?) defined by 

D{T*) = {u e L^{R^) I 3w„ e L^{R^) such that {vu\w) = {u\Tw) Vw e D{T)] 
T*u = Vu {vu is uniquely defined since D{T) is dense in L^(R'^)). 

The operator T is called self-adjoint if T* = T (i.e. if D{T*) = D{T) and if for all u e D{T) = 
D{T*), T*u = Tu). 

Let r be a self-adjoint operator on Li^iR?) with domain D(T), and z e C. In order to simplify 
the notation, we denote hj z — T the operator zlL2(jg3-) — T where 7^2(^3) is the identity operator 
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on L^(R^). If z — T is an invertible operator from D{T) to L^(]R'^), it can be proved that R{z) = 
{z — T)~^ defines a continuous linear operator on L'^{M.^) (with range D{T)). The set 



p{T) = {ze C I z-T is an invertible operator from D{T) to 

is called the resolvent set of T, and the family {R{z))z(^p{T) the resolvent of T . The spectrum 
of T is the set o{T) = C \ p{T). The set p{T) is an open set of C and o-(T) is a closed subset 
of K. An eigenvalue of T is a complex number A for which there exists u G such that 

Tu = Xu. The set of all the eigenvalues of T is called the point spectrum of T and is denoted by 
tTp(T). Obviously, <Jp{T) C cr{T) (in particular, all the eigenvalues of a self-adjoint operator are 
real). The set (Jc{T) = cr(T) \ crp(r) is called the continuous spectrum of T. If the continuous 
spectrum of T is empty (i.e. if a{T) — ap{T)), there exists a Hilbert basis (e„)„gN* of L^(R^) 
which diagonalizes T: 

T = ^ A„(e„| • )e„ = ^ A„ |e„) (e„| in bra-ket notation , 



with A,i e M. In this case (t{T) = {A,i}. If / : K — > C is continuous in a neighborhood of o-(r), 
the operator f{T) is defined as 

5^(l + |/(A„)n|u„P<oo| 



D(/(T)) = L=^M„e„e 2.2(^3) 

I riGN 



This definition can be generalized to any self-adjoint operator T by means of the Spectral Theo- 
rem [26j. Note that if T is non-negative, cf{T) C M+ and the operator T^/^ can therefore be given 
a sense. 

Lastly, the spectrum o'(T) of a self-adjoint operator can also be partitioned as follows 

a(r) = ad(r)U (T), 

where era (T) is the set of all the isolated eigenvalues of T of finite multiplicity, and where CTcss {T) = 
a{T) \ (Td(T'). The sets crd{T) and acss{T) are called respectively the discrete spectrum and the 
essential spectrum of T. 
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